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The transition between non resonant (Weibel-type) and resonant (whistler) instabilities is investi- 
gated numerically in plasma configurations with an ambient magnetic field of increasing amplitudes. 
The Vlasov-Maxwell system is solved in a configuration where the fields have three components but 
depend only on one coordinate and on time. The nonlinear evolution of these instabilities is shown 
to lead to the excitation of electromagnetic and electrostatic modes at the first few harmonics of the 
plasma frequency and, in the case of a large ambient magnetic field, to a long-wavelength, spatial 
modulation of the amplitude of the magnetic field generated by the whistler instability. 
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Electron distribution functions that are anisotropic in phase space are a common feature of collisionless plasmas 
both in space and in the laboratory and the investigation of the processes through which these distributions relax is 
of primary interest. In fact, the free energy that is made available by the unbalance of the particle "temperatures" 
in the different directions can be transferred, depending on the plasma conditions, to quasistatic magnetic fields, to 
electromagnetic or electrostatic coherent structures or to particle acceleration. 

The anisotropy of the electron distribution function in an unmagnetized plasma can give rise to the onset of the 
well known Weibel instability 1 which generates a quasistatic magnetic field. 

If a magnetic field is already present in the plasma, the Weibel instability driven by the anisotropy of the electron 
energy distribution turns into the so called whistler instability 2 in which case circularly polarized whistler waves 
are generated by the relaxation of the electron distribution function. Whistler waves are actually ubiquitous in 
plasmas and their generation has been extensively studied in recent years in the laboratory (see e.g. Ref|3|). Whistler 
instabilities have been reported in spaced where bursts of whistler mode magnetic noise are found to be present in the 
magnetosphere, close to the magnetopause and are also a likely source of several different magnetospheric fluctuations 
including plasmaspheric hiss and magnetospheric chorus. These waves propagate along the ambient magnetic field 
in the frequency range of Ud < uj r < w ce , where uj ce and oj C i are the electron and proton cyclotron frequencies. A 
I sufficiently large temperature anisotropy in a magnetized plasma is the most commonly observed mechanism for the 
l— "™ '■ generation of whistler waves, but additional mechanisms, such as trapped particle loss cone distributions, can also 
| generate these modes. 

■ The linear dispersion relations of these instabilities and how the Weibel instability merges into the whistler instability 
Q\ ' have been studied in the literature^. It has been shown^ that "Weibel-type" whistler modes occur in plasma conditions 
C*") . where the following independent inequalities are satisfied 
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where ui r and 7 are the real and imaginary parts on the mode frequency respectively, and depend on the parallel (v t h e ,\\) 
and on the perpendicular (vthe,±) thermal velocities of the electrons, and u> pe is the electron Langmuir frequency. On 
t-H . the contrary, if the above threshold, which we rewrite as 
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is passed, the growth rates of the whistler instability depart considerably from those of the Weibel instability. In this 
case the mode frequency uj r becomes comparable or greater than its growth rate 7. Here the anisotropy parameter A 
is defined by A = (v t he,±/v the ^) 2 - 1. 

In the present paper we consider the long term collisionless relaxation process of an anisotropic electron distribution 
focusing our attention on the nonlinear features of the transitional regime between the Weibel and the whistler 
instabilities and on the onset of secondary modes driven by the perturbations in the electron distribution function 
arising from the nonlinear development of the primary Weibel and whistler instabilities. 

The analysis presented here is performed by solving numerically the collisionless Vlasov equation for electrons and 
protons coupled to Maxwell's equations, in a restricted geometrical configuration (1D-3V) where all vector quantities 
are three dimensional but depend on one coordinate only, which is chosen to be the x coordinate, and on time t. 
More specifically, we consider a homogeneous initial plasma configuration with isotropic protons and a bi-Maxwellian 
electron distribution with T±_ > Tj | (where "parallel" refers to the x direction) in the presence of a uniform ambient 
magnetic field B x . The ratio between the proton temperature 7$ and the parallel electron temperature Tji is taken 
equal to one. Vector quantities must be taken to be three dimensional because whistler waves are circularly polarised 
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in the y-z plane and in addition we need to include longitudinal electric fields and velocities along x. This contrasts 
the treatment given for the Weibel instability (with ambient magnetic field B x = 0) given in Ref|6j where a f D-2V 
configuration was considered. 

At low ambient magnetic field, the nonlinear development of the whistler instability is similar to that of the Weibel 
instability, since the perturbed magnetic fluctuations are much larger than the ambient magnetic field. However, 
when the ambient magnetic field is sufficiently strong, the long term nonlinear behaviour of the instability changes 
drastically. In both limits the isotropization of the electron distribution function due to the onset of the instabilities 
is accompanied by the development of high frequency (Langmuir wave) electron density modulations^ 1 . These density 
modulations are forced by the spatial modulation of the magnetic energy density at wave numbers roughly two times 
the wave numbers of the most unstable Fourier components of the primary instabilities. However these density 
modulations become weaker as the ambient magnetic field is increased. In addition a purely kinetic effect is present 
in the small Bo limit where the development of the Weibel instability leads to strong deformations of the electron 
distribution function in phase space. These deformations have been shown in a 1D-2V configuration- to arise from the 
differential rotation in velocity space of the electrons around the magnetic field B z produced by the Weibel instability 
and to generate short wavelength Langmuir modes that form highly localized electrostatic structures corresponding 
to jumps of the electrostatic potential. These kinetic effects become weaker and eventually disappear as the ambient 
magnetic field is increased. 

In the case of the whistler instability, low-frequency density modulations involving the proton dynamics occur and 
generate soliton type structures known as whistler oscillitons^. Oscillitons are coherent nonlinear structures that occur 
in dispersive media for which the curves of the phase and group velocities cross for a finite value of the oscillation 
wavenumber k, see RefH Whistler oscillitons are of special importance as they have been invoked in order to describe 
coherent wave emission in the whistler frequency range observed in the Earth's plasma environment^. 



II. NUMERICAL SET UP 



In the numerical solutions of the Vlasov-Maxwell system presented below, all the parameters are normalized as 
follows. We use the plasma frequency uj pe and the velocity of light c as characteristic frequency and velocity; therefore 
the electron skin depth c/uj pe is the characteristic length scale. The electric and magnetic field are normalized to 
E = B = m e cLO pe /e and, finally, the electron temperature is normalized to m e c 2 . 

The initial electron anisotropy is kept fixed, T y = T z = Tj_, X] = T x , Tj_/Tj| = 12, while the value of B x is 
varied such that the corresponding electron cyclotron frequency, normalized on the electron Langmuir frequency and 
indicated in the following simply by Bq, changes from Bq = to Bq = 0.5. In these numerical integrations Ref|9| the 
ratio /3|| = v^ he ||/-E>o plays the role of a control parameter. Here the parallel thermal velocity v t he,\\ is normalized to 
c and is chosen to be equal to 0.02. The following values of B are chosen, B = 0, 0.02, 0.1, and 0.5, corresponding 
to /3n = oo, I, 4 x 10 -2 , and 1.6 x I0~ 3 respectively. 

The number of grid points in coordinate space is N x — 1024 and the length of the system is L x = IOtt in units of the 
electron skin depth d e = c/u pe . The total number of grid points in velocity space is N Vx — 101, N v ~ N Vz = 71. The 
system is evolved up to the time t = 1200 io~^ with time step 10~ 3 < dt < 10~ 2 . The total number of perturbed modes 
at t = is N = 220. The spatial grid spacing is dx — 0.03. Henceforth k max \]j <C 1 and k max d e ~ 1, where k max 
is the the most unstable wave number and Ad is the electron Debye length computed with the temperature Tji. The 
numerical integrations are initiated in a homogeneous plasma configuration with isotropic protons and Maxwellian 
electrons with different parallel and perpendicular temperatures by imposing a low amplitude magnetic field noise, 
randomly oriented in the y-z plane at different x positions: 
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B z (x,t = 0) = ^2 b ™ [cos{k n x + </>„)], (3) 

n=l 

with a,i n and bi n the initial "small" (~ 10 -4 ) amplitudes, ip n and <j) n random phases, n the wave number and 
k„ = 1ixnjL x the corresponding wave vector. 

In Table Q] we list the main parameters of the numerical runs, together with the analytical (kinetic) value, 7^™^ 
of the growth rate of the most unstable mode, k — k max , as obtained, after some simplifications, from Ref|5j, the 
corresponding value, 7„"™! as obtained by the numerical integration of the Vlasov-Maxwell system and the value 
of the ratio B z between the amplitudes of the z component of the magnetic field produced by the instability at 
saturation and that of the ambient magnetic field along x. In Table |TT] the real part of the frequency of the whistler 
unstable mode is shown for different values of k in the "large" ambient magnetic field case Bq = 0.5. For each k 
three values of the real part of the frequency are shown. Two are analytical: u>^ na (kinetic) is obtained, after some 
simplifications, from ReflU where a kinetic approach is used, and uif na (fluid) is obtained from Reffiol where a fluid 
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Table I: Main parameters of the numerical runs, the corresponding numerical and analytical growth rates for the most 
unstable mode corresponding to k = kmax and the normalized amplitude at saturation of the z component of the magnetic 

field produced by the instability. 



approach is used, while uj^ um is obtained by numerical integration of the Vlasov-Maxwell system. The comparison 
between the analytical and the numerical frequencies and growth rates shows good agreement. The numerical values 
of the frequencies and growth rates shown in Tables QUI] are obtained by considering only the linear part of the 
numerical integration. It should be noted that the comparison for the real frequencies is shown only for Bo = 0.5. 
For the other values of Bq in Table HI the oscillation period would turn out to be longer than the total duration time 
of the linear regime and it is thus difficult to estimate numerically. 
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Table II: Mode frequency as a function of k for Bo = 0.5 - Numerical and analytical values. 

In Table Hill the real part of the frequency of the most unstable whistler mode, as derived analytically in RefllOl is 
shown for different values of the ambient magnetic field. 
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Table III: Frequency of the most unstable whistler mode for different ambient magnetic fields 



III. WHISTLER MODES: TRANSITION FROM WEAK TO STRONG AMBIENT MAGNETIC FIELD 

From Table Q] we see that the growth rate of the whistler instability decreases as the amplitude of the ambient 
magnetic field increases, while the mode number k max , at which the maximum growth rate of the whistler instability 
is obtained, increases. 

By referring to the inequalities ([1} we see that for small ambient magnetic field amplitudes both inequalities are 
satisfied so that the mode frequency and its growth rate are approximately those obtained in RefJ^ for the Weibel 
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instability (zero ambient magnetic field). On the contrary, as the ambient magnetic field amplitude is increased, the 
system goes through the threshold condition ^ and at the same time the frequency of the whistler waves becomes 
comparable to, or larger than, their growth rate (uj r > J m ax)- Hence these numerical runs describe the transition 
regime between Weibel-like, low ambient magnetic field, instabilities and strong ambient magnetic field whistler 
instabilities. 



A. Electromagnetic fields 



The frequency spectrum of the z component of the perturbed magnetic field is shown in FigfT] for 
Bo = 0, 0.02,0.1, 0.5. For each value of the ambient magnetic field Bq the most unstable spatial Fourier mode 
(corresponding to k = k max in Table H| is chosen. The frequency spectrum is obtained by integrating over time over 
the whole simulation time (t = 1200 ^V, 1 ). We see a strong peak at uj r = 0.45 in the spectrum of the magnetic field 
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Figure 1: Frequency spectrum of the magnetic field component B z .k for Bq — 0,0.02,0.1,0.5 at the corresponding most 
unstable mode k = k ma x = 1.0, 1.2, 1.2, 2.0. Note the change in the amplitude scale in the different frames. 



with Bq = 0.5 (bottom right frame) which corresponds to the frequency of the whistler wave for k = k max — 2. This 
wave is observed at frequencies oj r = 0.015 and oj r = 0.07 for Bo = 0.02 and for Bq = 0.1, respectively. The whistler 
peak gets sharper, the stronger the ambient magnetic field. This reduction of the mode "line width" is consistent 
with the decrease of the mode growth rate at larger ambient magnetic field values. For comparison, the spectrum for 
the case B = (Weibel case) is also shown: no peak is observed as expected since in an unmagnctizcd plasma the 
instability is purely growing. We observe that both in the Weibel and in the whistler cases the frequency spectrum 
contains waves that propagate with a frequency equal to the plasma frequency, uj pe = 1 in the dimcnsionlcss units used. 
In the two intermediate cases, -Bo = 0.02 and Bo = 0.1 we also see waves at 2oj pe . Wave emission with characteristic 
frequency close to twice the electron plasma frequency are observed near the Earth's bow shocks and also from type 
II and III solar radio bursts 12 . Different mechanisms for the generation of such kind of wave emissions have been 
discussed in the literature, see Refill and references therein. Based on the frequency spectrum of the longitudinal 
electric field E x shown in FigO a plausible mechanism for the 2uj pe emission in our system can be the coupling be- 
tween forward and backward propagating Langmuir waves that generate an electrostatic mode at 2u; pe which is then 
converted into an electromagnetic wave by a nonlinear mode coupling process. The amplitude of the 2uj pe component 
of the perturbed magnetic field is much smaller than that of the waves radiated at the plasma frequency. This small 
amplitude electromagnetic emission at frequencies close to the harmonics of the plasma frequency is not limited to 
the modes at k = k max shown in FigQ] and extends over a range of different values of k. Within this range emission 
at 3uj pe is also observed (Figj3]) at short wavelengths. 

Electromagnetic emission at 3w pe has been reported in type II bursts^ and in type III bursts^. Various theories have 
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Figure 2: Frequency spectrum of the electrostatic field component E x> k for Bo = 0, 0.02, 0.1, 0.5 at the corresponding most 
unstable mode k = k maT — 1.6, 2.4, 2.8, 3.6. Note the change in the amplitude scale in the different frames. 
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Figure 3: Frequency spectrum of the Fourier component of the magnetic field, B z ,k for Bo = 0.1 at k = 1.4, 4.0, 5.0. Note the 

change in amplitude scale. 



been put forward in order to explain this higher-harmonic emission. The coalescence of a Langmuir wave and a 2w pe 
electromagnetic wave was suggested for the 3u; pe emission in RefH 

Furthermore, electromagnetic waves with a frequency u> r ~ k can be seen in the spectra in Fig|3](see also bottom right 
frame in FigJT]): their frequency does not correspond to harmonics of the plasma frequency and shifts with k. The 
amplitudes of these latter waves are smaller than those of the uj pe waves, but greater than those of the 2w pe and 3w pe 
emissions. The higher harmonics emission and the waves at uj r ~ k are present only at higher wave mode numbers. 
Multiple harmonic emission up to fifth harmonic in the Earth's fore-shock environment has been reported (see RefllTT). 
In our simulations the relative intensity of the emission decreases rapidly with increasing harmonic number with 
harmonics greater than three being too faint to be observed. Although in the magnetic field spectrum we do not 
observe electromagnetic emission at the fifth harmonic, such high harmonics are observed in the electrostatic spectrum 
as shown in Figf?l These uj pe and higher harmonic emissions are suppressed for larger values of the ambient magnetic 
field. In a typical solar wind configuration the ratio oj ce /u! pe is of the order— of 10 -2 . Hence only the choice of 
Bq = 0.02 or 0.1 represents a fairly typical situation for this kind of emission. This can be clearly seen from Fig[S] 
where a contour plot of the perturbed magnetic field B z in the uj r -k plane is shown. 
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Figure 4: Frequency spectrum of the Fourier component of the electrostatic field, E x ^ at k = 1.4,4.0,5.0 for Bo = (top 

frames) and Bo = 0.1 (bottom frames). 



B. Electrostatic fields 



The frequency spectra of the longitudinal field E Xt k at k — 1.4, 4, 5 are shown in Fig[3]for Bq = and for Bq = 0.1. 
In both cases the spectra exhibit a clear peak at the electron plasma frequency. This peak for Bq — appears to be 
less pronounced for larger values of k, consistently with what shown in Ref|6|, while for Bq = 0.1 the amplitude at 
tdpe remains almost the same even at smaller wavelengths. In addition in both cases we also see peaks appearing at 
multiple harmonics of the plasma frequency. It should be noted though that for the unmagnetized plasma case only 
harmonics up to 4w pe are observed in contrast to harmonics up to 5uj pe that are observed in the case with Bq = 0.1. 
The amplitudes of the multiple harmonics of for Bq = 0.1 are much larger than in the B = case. In addition, 
similarly to the 1D-2V, Bq = analyzed in Reflg, electrostatic structures are formed corresponding to plasma modes 
with values of k larger than those of the whistler waves. For example, for Bq — 0, the resonant velocity corresponding 
to the clear bump in the distribution function shown in Fig|8](top right frame) is v r = 0.05, leading to the excitation of 
Langmuir waves with wave mode k ~ uj r /v r ~ 20, while for Bq = 0.1, the resonant mode at which the Langmuir waves 
are excited is k ~ 10. As discussed in Ref@, these modes are excited by the deformation of the electron distribution 
function due to the differential rotation of the electrons in phase space in the magnetic field generated in the y-z 
plane by the Weibel and by the whistler instabilities. This deformation becomes weaker as the effect of the ambient 
magnetic field Bq is increased. In addition if we compare the results obtained for Bq = in the 1D-2V case (see Fig. 12 
in Ref|(|) with the corresponding results shown in Fig(S] (left frame) for the 1D-3V case, we see that in the latter case 
the amplitude of these structures is somewhat smaller than in the 2V case. This difference can be attributed to the 
fact that the 2V case imposes a higher coherency of the perturbed magnetic field which is necessarily oriented along 
the z axis. On the contrary, in the 3V case considered here, the orientation of the perturbed magnetic field in the y-z 
plane can change as a function of x. The eventual formation of coherent electrostatic structures can be attributed 
to the fact that, for sufficiently long times, the initially randomly orientated perturbed magnetic field becomes fairly 
organized, as shown by FigjB] (right frame) where the angle a — arctanh(B y /B z ) is shown as a function of space and 
time. This allows for the formation of peaks in the electron distribution function, that are sufficiently coherent in 
space so as to excite plasma waves by inverse Landau damping, see Figs 13 and |S] The contour plot of the electron 
distribution function in the v x -v y plane for Bq — is shown in FigJTl The distribution function starts being distorted, 
with a differential rotation combined with a spreading along v x , at the time (t — 300) when the instability begins to 
saturate. These deformations become 'multi-armed' as time evolves. Similar kind of deformations are observed in the 
v x -v z plane (not shown here). Multi-armed structures are also noticeable in the v y -v z plane (FigJTl bottom frames). 
The combined velocity and space dependence of the electron distribution function is shown in Fig[8l bottom frames. 
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Figure 6: Left frame: contour plot of the electrostatic field in the x-t plane (note only a part of the time interval of the total 
time of the simulation is shown). Right frame: the angle a = arctanh(B y / B z ) between the magnetic field components in y 

and z direction in the x-t plane. 



The left frame shows the formation of phase space vortices in the x-v x plane (shown at v y = and v z = 0). The 
right frame shows alternated filamented structures in the x-v z plane. The distribution function for Bo = 0.1 behaves 
similarly (not shown here), but the strength of the deformations decreases. 



C. Low frequency modulations 

As we have already discussed in Sec lIIIBl the distorted distribution function for Bq — 0.1 leads to the resonant 
excitation of short wavelength plasma waves at large wave numbers, e.g. (fc ~ 10). On the contrary, at smaller 
wave numbers we see the interplay between the inhomogeneity of the perturbed magnetic field amplitude and of 
the electron and proton densities. The electrostatic and the magnetic fields are shown in x-t space in Fig[5] for 
Bo = 0.5. These figures show the forward and the backward propagation of the whistler waves. In addition, long 
wavelength modulations of the electrostatic structures and of the proton and electron densities are seen in FigfTU] at 
half the wavelength of the perturbed magnetic field. The wavelength of the modulations of both the magnetic field 
and of the proton and electron densities decreases with increasing ambient magnetic field, the electrostatic field and 
density modulation wavelength still being equal to half the wavelength of the most unstable Fourier component of the 
perturbed magnetic field. We focus our attention primarily on the case with ambient magnetic field Bq = 0.5. 
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Figure 7: Top frames: Contour plot of the electron distribution function in the v x -v y plane (for v z = 0) for Bq = at 
x — 1.86 (position corresponding to the maximum value of the z component of the perturbed magnetic field) for t = 250, 
t — 300, t = 350. Bottom frames: Contour plot of the distribution function in the v y -v z plane (for v x = 0) at x — 1.86. Note 

the different velocity ranges. 
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Figure 8: Top frames: The normalized electron distribution function versus v x at x = 1.86 for v y = and v z — for Bo = 0. 
Bottom (left) frame: Contour plot of the distribution function in the x-v x plane at v y — and v z — and x = 1.86. Bottom 
(right) frame: Contour plot of the distribution function in the x-v z plane at v x = and v y = and x = 1.86. 



There are two kinds of modulations in the perturbed magnetic field, one with long wavelengths having a width almost 




half the total length of the simulation box, L ~ 17 and corresponding time period T ~ 200. The phase velocity of 
these modulations is v p h ~ 0.08. The other kind has shorter wavelengths, i ~ 2, and T ~ 7, i.e. u p h ~ 0.28. The 
range of phase velocities corresponding to the peak in the whistler waves frequency spectrum for different values of k 
extends from 0.005 < v p h to v p h < 0.25. Hence, the long wavelength phase velocity u p /j ~ 0.08, as obtained from the 
magnetic field modulations, falls in the range of the whistler phase velocities and it can be further noticed that it is 
comparable with some of the phase velocities obtained for both short wavelength and long wavelength electrostatic 
modulations which are in the range 0.005 < v P h < 0.01 and 0.075 < v p h < 0.13, respectively. Furthermore the 
range of phase velocities 0.08 < v p h < 0.01 obtained from low frequency peak, 0.03 < u) r < 0.04 in the frequency 
spectrum of the electrostatic field at the most unstable mode k max = 3.6 (see Figf2] bottom right frame) covers a 
large part of the phase velocity range of the whistler waves. We can thus infer that the dominant whistler waves 




Figure 10: Electron density and proton density for Bo = 0.5. 



are modulated by the low frequencyperturbations observed in the electrostatic field. This kind of modulations has 
been studied theoretically. In Reffl9l it was proposed that the nonlinear coupling of oppositely travelling waves 
(forward and backward waves) produce a low frequency ponderomotive force with frequency u — lu\ — ui2 and wave 
vector k = k\ — k-2 = 2k\. In Refs. [2^ and l2ll it was later shown that these ponderomotive forces are reinforced 
due to the interaction of the whistler packets with the background low density perturbations — . In addition to the 
low density perturbations observed in our simulations (see Fig [10]), the formation of wave packets of the perturbed 
magnetic field are also evident (see Fig llll left frame). Coherent wave emission in the whistler frequency range 
(0.1 — 0.5w ce ), consisting of nearly monochromatic wave packets, has recently been observed by several spacecraft in 
the Earth's plasma environment^. These are mainly generated by oscillitons which are stationary, nonlinear structures 
exhibiting solitary structures with embedded smaller-scale oscillations resembling wave packets. These arise from the 
momentum exchange between protons and electrons. These oscillitons in the whistler frequency range are known as 
'whistler oscillitons'. The appearance of stationary nonlinear waves is related to the existence of a 'resonance point' 
where phase and group velocity coincide. The right frame in Fig llll shows the range of frequencies and wave number 
available from the numerical results which have a value of their phase velocity close to that of the corresponding group 
velocity. We also observe a shift of the whistler waves to longer wavelength corresponding to the shift of the broad 
magnetic spectrum to lower k values during the saturation regime when the coherent whistler emission takes place 
and modifies the frequency range of the whistler modes. 
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Observations from laboratory experiments^ exhibit a clear evidence of this kind of modulated whistler wave packets 
due to nonlinear effects. Furthermore, instruments on board CLUSTER spacecraft have observed^ broadband intense 
electromagnetic waves. Cluster measurements also exhibit the formation of envelop solitary waves accompanied by 
plasma density cavities. Both electron and proton density cavities formed in our numerical results are shown in Fig llOl 
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Figure 11: Left frame: Magnetic field versus position at t = 540. Right frame: The dispersion relation of the whistler waves 
available which have a value of their phase velocity close to that of the corresponding group velocity. 



IV. CONCLUSIONS 



The anisotropy of the electron distribution function ( Tj_ > Xj|) has been shown to lead to the onset of both resonant, 
w r > "/ m ax, and non resonant, io r « 7 ma£ c, instabilities depending on the value of the ambient magnetic field, see 
Table HI When the plasma is unmagnetized the instability (Weibel instability) is completely nonresonant, where 
uj r = << j max — 0.038 (in units of the Langmuir frequency), but as we move towards higher ambient magnetic 
fields, ui r becomes comparable to the maximum growth rate of the instability. At very high ambient magnetic field 
(Bo = 0.5, being the value of the electron cyclotron frequency normalized to the Langmuir frequency) the mode 
becomes resonant with uj r much larger than the growth rate (whistler instability). 

The development of these instabilities is accompanied by the generation of electrostatic and electromagnetic waves at 
the harmonics of the plasma frequency (up to the fifth harmonic in the case of the electrostatic spectrum). 
In addition the spatial self-organization of the perturbed magnetic field leads to the deformation of the electron 
distribution function and to the consequent generation of short wavelength electrostatic modes. These secondary 
phenomena tend to be suppressed as the ambient magnetic field becomes larger. On the contrary, at large ambient 
magnetic field we find a low frequency long wavelength modulation of the whistler wave spectrum that we interpret 
in terms of the self-consistent interaction between the perturbed magnetic field and the low frequency electron and 
proton density modulations. The magnetic modulations that we observe are reminiscent of the so called "whistler 
oscillitons" that arise under the condition that the mode phase velocity coincides with its group velocity. 

The simulation results presented here provide relevant information on the different processes that the presence of 
anisotropy can drive in collisionless plasmas in different plasma magnetization regimes, as parametrized by the value 
of the parameter j3±_ > 1 that gives the ratio between the perpendicular plasma pressure and the ambient magnetic 
field pressure. For example, in the context of solar physics, the occurrence of resonant whistler modes has been 
advocated as an effective acceleration mechanisms of electrons directly from the background plasma up to ~ 25 MeV 
within several seconds under normal solar flare conditions^. 

The simulation results presented in this paper have been obtained by using a ID configuration and are thus restricted 
to the study of the nonlinear effects on the system only of modes with parallel wavenumber k. This restriction will be 
lifted in a future paper. However, linear analysis^ indicates that in the presence of low ambient magnetic fields the 
whistler instability remains the fastest growing one even in a two dimensional configuration where we allow modes to 
propagate at an angle with respect to the ambient field. On the contrary, at higher ambient magnetic fields where 
the whistler instability is suppressed, obliquely propagating modes have a larger growth rate. An important question 
to address in this latter case is how the formation of electrostatic structures will develop in such a regime. 
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